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Abstract 

When a porous rock is saturated with an electrolyte, electrical fields are 
coupled with seismic waves via the electro-seismic conversion. Pride [6] derived 
the governing models, in which Maxwell equations are coupled with Biot's 
equations through the electro-kinetic mobility parameter. The inverse prob- 
lem of the linearized electro-seismic conversion consists in two step, namely 
the inversion of Biot's equations and the inversion of Maxwell equations. We 
analyze the reconstruction of conductivity and electro-kinetic mobility pa- 
rameter in Maxwell equations with internal measurements, while the internal 
measurements are provided by the results of the inversion of Biot's equations. 
We show that knowledge of two internal data based on well-chosen boundary 
conditions uniquely determine these two parameters. Moreover, a Lipschitz 
type stability is proved based on the same sets of well-chosen boundary con- 
ditions. 



1 Introduction 

When a porous rock is saturated with an electrolyte, an electric double layer is 
formed at the interface of the solid and the fluid. One side of the interface is 
negatively charged and the other side is positively charged. Such electric double 
layer(EDL) system is also called Debye layer. Due to the EDL system, electro- 
magnetic(EM) fields and mechanical waves are coupled through the phenomenon 
of electro-kinetics. Precisely, electrical fields or EM waves acting on the EDL will 
move the charges, creating relative movement of fluid and solid. This is called 
electro-seismic conversion. Conversely, mechanical waves moving fluid and solid 
will generate EM fields. This is called seismo-electric conversion. Thompson and 
Gist [S] have made field measurement clearly demonstrating seismo-electric conver- 
sion in saturated sediments. Zhu et al. [TTl [TH [TJ] made laboratory experiments 
and observed the seismo-electric conversion in model wells, and their experimen- 
tal results confirm that seismo-electric logging could be a new bore-hole logging 
technique. 

The investigation of wave propagation in fiuid-saturated porous media was early 
developed by Biot [31 [3] . The governing equations of the electro-seismic converstion 
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was derived by Pride fS' as following. 



V X E ^ iuj^H, 
V X = (cr - ieoj)E + L{-Vp + oj'^pfu) + J, 

— UJ^{pU + PfW) = V • T, 

-iLuw = LE + ^{-Vp + uj'^pfu), 
T = (AV • u + cV • w)I + G{Vu + Vu^), 
-p = cV • M + MV • w, 



s J 



(1) 

(2) 
(3) 
(4) 

(5) 
(6) 



where the first two are Maxwell's equations, the remaining are Biot's equations. 
The notation is as follows: 



E 


electric field, 


H 


magnetizing field or magnetic field intensity. 




seismic wave frequency. 


a 


conductivity, 


e 


dielectric constant or relative permittivity. 


A* 


magnetic permeability. 


Js 


source current. 


P 


pore pressure. 


Pf 


density of pore fluid. 


L 


electro-kinetic mobility parameter. 


K 


fluid flow permeability. 


U 


solid displacement. 


W 


fluid displacement. 


T 


bulk stress tensor, 


V 


viscosity of pore fluid. 


X,G 


Lame parameters of elasticity. 


C,M 


Biot moduli parameters. 



Pride and Haartsen [7] also analyzed the basic properties of seismo-electric waves. 

Notice that the coupling is non-linear, namely electro-seismic and seismo-electric 
conversions happen simultaneously. Under the assumptions that the coupling is so 
weak that multiple coupling is neglectable, we can linearize the forward system 
in two steps. Particularly, we focus on the electro-seismic conversion and ignore 
the seismo-electric conversion. The first step in the forward system is modeled 
by Maxwell equations without the effect of the seismic waves, i.e., L = in 0. 
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While the electro-seismic conversion happens, the seismic waves are generated and 
modeled by Biot's equations with potential LE in (jj]). 

In the present paper, we mainly focus on the inverse problem of the linearized 
electro-seismic conversion, which is a hybrid problem and consists of two steps. 
The first step of the inverse problem is to invert Biot's equations, i.e., to recover 
the potential LE in (jH) from any measurements observed on the domain boundary. 
Williams [10] presented an approximation to Biot's equations, which could be a 
useful tool to study the inverse problem. 

Assuming the first step is implemented successfully, the second step of the in- 
verse problem is to invert Maxwell's equations, which consists of reconstructing the 
conductivity a and the electro-kinetic mobility parameter or the coupling coefficient 
L from boundary measurements of the electrical fields and the internal data LE 
obtained in the first step. 

The problem of interest in this paper is the second step of the inverse problem. 
We study the reconstruction of the conductivity a and the coupling coefficient 
L and prove uniqueness and stability results of the reconstructions. Particularly, 
we show that cr, L are uniquely determined by 2 well-chosen electrical fields at the 
domain boundary. The explicit reconstruction procedure is presented. The stability 
of the reconstruction is established from either 2 measurements under geometrical 
conditions or from 6 well-chosen boundary conditions. 

Mathematically, our proof relies on explicit solutions to Maxwell's equations, 
namely Complex Geometrical Optics (CGO) solutions, constructed by Colton and 
Paivarinta[5]. In our reconstruction procedure, the coupling coefficient L satisfies 
a transport equation with vector field /3. With CGO solutions, we can prove the 
integral curves of the vector field /? are close to straight lines and exit the domain in 
finite time. Therefore, L can be uniquely and explicitly solved by the characteristic 
method. Stability follows the analysis of the method of characteristic. 

The rest of the paper is structured as follows. Section[2]presents our main results. 
The CGO solutions are introduced in section [01 The inverse Maxwell's equations 
and an explicit reconstruction algorithm are addressed in the rest of section [31 while 
section 13.21 focusing on the proof of the uniqueness result and section 13.31 and 13.41 
focusing on the stability proof. 



2 Main results 

Let ft be an open, bounded and connected domain in M'^ with boundary dfl. In 
the second step of the electro-seismic conversion, the propagation of the electrical 
fields is modeled by Maxwell's equations in il, 

V X E = iujfiH, , , 

V X = {a -ieiu)E + Js. ^ ' 

The measurements available for the inverse problem include the internal data from 
the first step 

D := LE, in ^ (8) 

and the boundary illumination, i.e., the tangential boundary measurement of the 
electrical field 

G := tE, on (9) 
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Define the operator 

AAf(L,a) := (J^AG). (10) 

The problem now is to invert the operator Am, or namely, to reconstruct {L,a) 
from some measurements {Js,j,Dj, Gj) indexed by j, assuming ^, and e are given. 

The main purpose of this paper is to prove the uniqueness and stability of the 
coefficient reconstructions. Wc define the set of coefficients (L, cr) € as 

M = {{L, ajeC X C^-^ ■.d>2 (11) 
and is not an eigenvalue of V x V x • — fc^n}, 

where the wave number fc > and the refractive index n are given by 

k = ioWeofio, n= — fe + i-V (12) 

The main results are as follows, where the measurements G and D are complex- 
valued. 

Theorem 2.1. Let H, be an open, bounded subset of with boundary dfl of 
class C^. Let {L,a) and {L,a) be two elements in M. Let D := {Di,D2) and 
D := (Z)i,I?2), be two sets of internal data on for the coefficients {L,a), {L,a), 
respectively and with boundary illuminations G := (Gi,G2). 

Then there is a subset of G & (C"'+^(9ri))^, such that if Dj = D^, j = 1,2, we 
have {L,a) = {L,a). 

To consider the stability of the reconstruction, we need to restrict to a subset of 
f2. Let Co be a constant unit vector. Let .xq G dVl be the tangent point of dVt with 
respect to Co> i-e-, the tangent line of dVL at Xq is parallel to Co- Define i\ to be the 
subset of f2 by removing a neighborhood of each tangent point xq e 90. 

Theorem 2.2. let d>2>. Let Q. be convex with C"^ boundary dVl and Oi is defined 
as above. Assume that {L,a) and {L,a) are two elements in A4. Let D = (Dj) 
and D = (Dj), j = 1,2, be the internal data for coefficients (L,a) and {L,a), 
respectively, with boundary conditions G — {Gj), j — 1,2. 

Then there is a set of illuminations G e (C"'+^((?0))^ such that restricting to 
Oi, we have 

\\L - L||c"'-i(Oi) + Ik - o-||cd-3(f2i) < C\\D - D\\(C''(niW- (13) 

The geometric conditions can be removed when more measurements are avail- 
able. In particular, when 6 complex measurements are provided, we have the fol- 
lowing stability result. 

Theorem 2.3. let d > 3. Let n be convex and Qi is defined as above. Assume that 
{L,a) and {L,d) are two elements in Ai. Let D = {Dl,D^) and D = {D{,D2), 
j = 1,2,3, be the internal data for coefficients {L,a) and {L,a), respectively, with 
boundary conditions G = {G{,G2), j = 1,2,3. 

Then there is a set of illuminations G € (C'^'^^ (dil))^ , such that 

\\L - LWc^-iin) + [\a - a\\c^-3^n) < C\\D - -D||(cd(o))6. (14) 

Note that the above measurements are all complex-valued. We will need two 
real measurements to make up one complex data. 
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3 Inversion of Maxwell's equations with internal 
data 

Let rt be an open, bounded and connected domain in M.^ with boundary dfl. In 
the case when /i = /iq is constant and = in Q, we can rewrite the system in ([7]) 
as 

V xV X E - k'^nE ^ 0, (15) 

and 

V ■nE = (16) 
where the wave number fc > and the refractive index n are given by (jl2p . 



3.1 Complex Geometrical Optics solutions 

Colton and Paivarinta_5j constructed exphcit solutions, namely Complex Geomet- 
rical Optics solutions(CGOs), to the Maxwell's equation (|T5|) and (|T6)) . CGOs will 
be the main technique we will use to solve the inverse Maxwell problem. We follow 
the construction of CGOs in [S] and extend their properties from space to higher 
order sobolev spaces. CGOs are of the form 

Eix) = e^i-^i7^ + Rcix)), (17) 

where C £ C''^\M'', 77 G C'^, are constant vectors satisfying 

C•C = A;^ C-7? = 0. (18) 

Substituting ^ into ^ and ^ gives 

V X V X = k^{n-l)-q + k^nRc^, (19) 
V-i?c ^ ~a-{T] + Rc) (20) 

where V := V + iC a :— Vn{x) /n(x). We further define A := A + 2i( • V — fc^. 
By substituting the formula V x V x i?^ = -ARq + VV • Rq into HI]) and (HOI), we 
see that i?^ is a solution to 

(A + 2< • V)i?c = -V(a • (r; + R^)) + k^l - n){ri + R^;). (21) 

It was proved in [S] the existence of i?^ to pT|) as a C^(M'^) functions. For our 
analysis, we need to extend the results of CGOs in [5] to smoother function spaces. 

Let the space Lg for (5 G M be the completion of C^(R^) with respect to the 
norm || • ||^2 defined as 

\\u\\^. = (^ljxr\u\'dx^ \ {x) = {l + \x\^)i. (22) 

To get smoother CGOs than that constructed in 5^, we introduce the space Hg for 
s > as the completion of C^(M'^) with respect to the norm || • ||^f| defined as 

\\u\\Hi = (^Jjxf'\{I-A)iufdxy (23) 
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Here (/ — A)iu is defined as the inverse Fourier transform of where ■u{^) 

is the Fourier transform of u{x). 

We recall [S] for |CI > c > and v e Lj^^ with — 1 < (5 < 0, the equation 

(A + 2< ■y)u = v (24) 

admits a unique week solution u £ L'j with 

\\u\\^.<C{6,c)^^. (25) 

Since (A+2iC-V) and (/—A)" are constant coefficient operators and hence commute, 
we deduce that when i; g for s > 0, then 

lklk|<C(5,c)^^. (26) 



We define the integral operator G(; : HI^^{W) HI{W) by 



where is the inverse Fourier transform. We see that is bounded and there 
exists a positive constant C{6) such that 

IIGcll < ^. (28) 

Before we can prove the existence of a unique solution to (j2ip . we first prove the 
following lemma. 

Lemma 3.1. For any v G -ff|_|_j^(K'^) and |C| sufficiently large, the equation 

(A + 2< • V + a • V)u = V (29) 
has a unique solution u G H^{^^) satisfying 

\\u + n-^'^G<,{n''\)\\Hi<^, (30) 

Jor some positive constant C independent of Q. 

The Lemma 3.1 in 5 proves for the case when s — Q. We study any s > here. 
Proof. From the identity 

n-i/2(^ + 2iC ■ V){n^l'^u) = (A + 2iC • V + a • V)u + qu, (31) 
where q := An^/^/n^/^, we can rewrite (|29|) as 

(A + 2zC-V-(7)/ = 5, (32) 
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where / := n}/'^u and g := n^l'^v. Applying the integral operator Gq gives 

/ + Gc(g/) = -Gc(g), (33) 

which admits a unique solution in i?|(]R^) since / + Gi^{q-) is invertible for \(\ 
sufficiently large. Eq. pS]) also gives 



11/ + Gd9)\\Hi = \\Gdq{Gdqf) + Gc(5)))IU| < (34) 

for some positive constant C independent of This proves the lemma. □ 

Proposition 3.2. For |C| sufficiently large, there is a unique solution G Hg{M.^) 
to (dl]). Thus, the COO solution E defined by dTT]) satisfies ([13 and (dH). More- 
over, R(^ satisfies 

\\R-^n-'/^Gcin'/'a■rJ)C\\H|=0(^■^^. (35) 

Proof. By applying the vector identity 

V{A ■ B) ^ Ax {V X B) + B X {V X A) + (A- V)B + {B ■ V)A, (36) 

and rearranging terms, we can rewrite (|2ip as 

(A + 2< • V + a • V)Rc = - a x {\/ x R^) - {R^^ ■ V)a (37) 

- V{a ■ r/) + fc2(l - n){T] + R^). 

Recall that V = V + iC, which is potentially troublesome. We assume that Q :— 
V X i?f . By dm), we have that 

y xQ^k^{n-l)T] + k^nRc (38) 

and hence 

V X V X Q = fc^Vn X (r/ + R,^) + k^{n - 1)( x r] + k^uQ. (39) 

Since V • Q = 0, we now have 

AQ 2< • Vg = /c^Vto x (yy + i?^) + ^^(l - n)« x r/ + Q). (40) 

If Ri^ £ Hg and \(\ is large, ([25)) implies that / + k^G(;{l — n) is invertible and the 
IIQIIi/l is bounded by a constant independent of C- 

Therefore, by applying Lemma lB.ll to p7p and substituting the solution Rq in to 
the right hand side of ([37]) recursively, Rq is given as a summation of a power series 
of Gij, which converges due to ([28]). This proves that (|37l) has a unique solution in 
H^{S?). Furthermore, by Lemma [3TT1 we see that 

\\R + n-^^'^Gc^(n^l\~a x (V x i?^) - (i?^; • V)a (41) 

- V(a • ^) + fc2(l - n){^ + i?c)])IU| = o (^^) . 
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Substituting V(a • 77) = (V + iC){a ■ 77), (gT]) implies that 

\\R+in-^^^G^in^/^a ■ r,)C\\Hi < \\n-^/''G^in^^^[~a x (V x i?^) (42) 
- • V)a - Va • 77) + fc2(i _ n)(r; + Rc)])\\hi + O ( 



,ICI 

This complete the proof. 



\C\' 

O ( ir) (43) 



□ 



By applying Sobolev embedding theorem to on a bounded domain, we have the 
estimate in _ff*(il) and C^, when s ^ ^ + d + l, l > 0. 

Corollary 3.3. Let be an open, bounded domain. With same hypotheses as the 
previous proposition, we then have 

- in-'/^Gc{n'/'a ■ v)C\\H^in) < j^, (44) 

for some positive constant C independent of C,. Moreover, when s = | + d + t, we 
also have 

||i?C - *n-i/2G^(„i/2c, . r;)C||c<<(o) < ||, (45) 
for some positive constant C . 

Proposition 3.4. Suppose C, G C'^\M'^, 77 G C^, satisfy ( ■ ( = k'^ and C • 77 = such 
that as I CI —>■ 00 the limits C/\C\ OLnd rj exist and, 

IC/IC|-Co|=o(^), \rj-r^o\ = o(^l-^. (46) 

E Hg{M.'^) is the unique solution to (pij) in Provosition VJ.SX Let s — | + (i + i. 
For ICI large, 

- i\C\n-^/^G^{n'/^a ■ iiMc-i^n) = O (^^) . (47) 

Proof. The proof follows directly by substituting into (H5|) . □ 

We choose the specific sets of C, 77 as in [5,. Precisely, Let ft, be a small real 
parameter and choose arbitrary a G R. We define C,i , C2 and 771 , 772 by 



Ci = (a/2,i^l/ft2 + a2/4-fc^f//7). 



C2 = (a/2, -i^l/h^ + a2/4 - P, -l/h), 

m = ^7^(i//7,o,-a/2), 

m = ^^(l/ft,0,a/2). 
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and note that 



limry, -r?o := (1,0,0), J = 1,2, 

c— >-oo 

^limCi/ICil = Co := 7|(0,«,1), (49) 

liniC2/|C2| = -Co, 



and 



Ci+C2 = (a,0,0), Co-Co = 0, 77o-Co = 0. (50) 
Proposition 13.41 implies that 

(7yi+i?cJ-(ry2 + i?c.)-l + o(l) (51) 
in the norm over bounded domain f2. 

3.2 Construction of vector fields and uniqueness result 

Let us now consider the reconstruction of {L,a). Assume Ej, j — 1,2, be two 
complex solutions to 

V X V X - k'^nEj = in f7, (52) 
with the tangential boundary conditions 

tEj = Gj on dVL, (53) 

with Gj well-chosen boundary values and j = 1,2. We will see that 

V X V X £^1 • - V X V X £^2 • -Bi = 0. (54) 

Let Dj = LEj, j = 1,2, be the internal complex-valued measurements. Assume 
L e C"'+^(il) is non- vanishing. Substituting Ej = Dj/L in we have, after 

some algebraical calculation, 

/3-VL + 7L = 0, (55) 

where 

P = xix)[{yDi)D2-{VD2)Di] + [(y -0^)02 ~{V -02)01] (56) 

-2[{S/DifD2~{WD2)0,], 
7 = x{x)[yi^ -Di) -02- ■O2) -01] + [S/^Oi -02-^/^02 -Oi]. {57) 

Here, x(a;) is a smooth known complex- valued function with |x(2;)| uniformly bounded 
from below by a positive constant on f2. 

To show the transport equation (|55p has a unique solution, it suffices to prove 
that the direction of the vector field /3 is close to a constant and thus the integral 
curves of /3 connects every internal point to two boundary points. 

Let El, E2 be two CGOs with parameters Ci, C2 and rji, 772 defined in (|48)) . i.e.. 



E,^e'iy-irj,+R^.), j = 1,2. (58) 
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Let Dj = LEj, j — 1,2, be the correspondmg internal data. By choosing x{x) = 
_Q-iiCi+C2)-x_^ and substituting Dj into (|56)) . we can analyze the asymptotic 

behavior of the vector field /3 as |Cj | — > oo, or equivalently, /i — >■ 0. Indeed, we have 
V^, = e<^--[(77, + Rc,){WLf + tfiir^, + i?c,)Cj + ^Vfe + R^)], J = 1, 2, (59) 
and 

[(7?i + i?cJ(VL)^ + iL(7?i + i?cJCf + MV(r/i + i?cJ](?72 + Rq) 



4\/2 

K'?! + i?cJ[VL^(r;2 + i?c J + ^^Cf (^/2 + i?cJ] (60) 



4^/2 

+mV(?7i+%)(?72 + %)] 

^im + Rci)[':i-iv2 + Rc.)] + 0{h). 



4\/2 

Therefore, by Proposition O x(a;)(V£»i)I)2 ^ in C"'(r2) norm as — >■ on a 
bounded domain f2. Similarly, 

x(a;)(VDi)^752 

= -^[VL(77i + i?ci)^(^2 + i?cJ + «iCi('?i + i?o)^('?2 + i?cJ 

+LW{rn + Rc,f{m + Rc2)] (6i) 

_^ inC^in) as /i^O. 

4 

More calculation gives that 

= -^[(V/i- (7?i + i?cJ)(r72 + i?c2) + ii(Ci • im + Rci))im + Rc.) (62) 

+LV-(r,i+i?cJ(^2+i?C2)] 
^ in C"^(f]) as h~^0. 

By substituting dBHl), (EB) and (HI]) into dSS]), we have 

lim||/3-i2^ol!c<'(o) =0, (63) 

i.e., the vector fields have approximately constant directions for small h and their 
integral curves connect every internal point to two boundary points. Thus, the 
transport equation (|55p admits a unique solution. 

To see the dependence of vector fields on the boundary conditions, we need to 
introduce a regularity theorem of Maxwell's equations. Let tE be the tangential 
boundary condition of E. Define the Div-spaces as 

H^i^{n) ^{ue H'n^il) : Div(tu) G H'-^/^{dn)}, (64) 
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THf^i^idn) = {g e H'n'^{dn) -. Div(g) e H'{dn)}, (65) 

where H^il.^{n) is a space of vector functions of which each component is in H'^(n). 
These are Hilbert spaces with norms 

\\u\\Hi„^{n) = + l|Div(tu)||^.-i/2(an), (66) 

\\g\\THf^.^{dn) = ||.g||H=(ao) + ||Div(.g)| 1^,(00) (67) 

It is clear that t{Hl).^{n)) = TH'^^^'^{dn). 

Theorem 3.5. Let e,/i e C", s > 2, he positive functions. There is a discrete 
subset S C C such that if oj is outside this set, then one has a unique solution 
E e to fTS]) given any tangential boundary condition G G TH^-^^^^dQ). The 

solution satisfies 

< C||G||^^.-v.(,o) (68) 

with C independent of G. 

Note that when the tangential boundary condition is prescribed by CGOs, i.e., 
Gj ^ tEj, j = 1,2. By Theorem [331 Ej is the unique solution to ([SI]) and ([5^ . 
Then the corresponding vector field $ defined in ([55]) satisfies which implies 
that the direction of /? is close to constant direction and thus its integral curves 
connect every internal point to two boundary points. Therefore, (|55p admits a 
unique solution. 

Furthermore, Theorem l3 . 51 also allows one to relax the boundary condition Gj = 
tEj and still to get the uniqueness of the solution to (|55|) . 

Proposition 3.6. Under the assumption of Theorem \3.5[ when Gj is in a neigh- 
borhood of G = tE in C'^^^{dn), j ~ 1,2, the corresponding vector field (3 defined 
in ([55]) satisfies 

l|/3-i'Co||c-(o) =0(/i), (69) 

for small h. 

Proof Let s = f + d + i. Notice that C^+^idO.) C H'^+^{dVt) C Hi+'^+'fd^). By 
the Sobolev embedding theorem and Theorem 13.51 we have that 

11-^11^+1(0) < C'l|G||cd+3(en)- (70) 

Let us now define boundary conditions Gj G C'^^^{dil), j — 1,2, such that 

ll'^j" ~ tE-j\\cd+:ii^Q^) < e, (71) 

for some e > sufficiently small, let Ej be the solution to (|15p and p6)) with 
tEj Gj. By jTO]), we thus have 

ll^i ~ -E^j 11^+1(12) < Ce, (72) 

for some positive constant G. Define the complex valued internal data Dj — LEj. 
We deduce that 

<Coe, (73) 

for Co > 0. Define P by We can easily deduce (jMl) from ([63]) and ([73|). This 

finishes the proof. □ 
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Recall Ai is the parameter space of (L, a) defined in ([TT|) and h is the parameter 
in (|48)) . We are in the place to prove Theorem 12. II 

Proof of Theorem \2.1[ By Proposition 13.61 we choose the set of illuminations as a 
neighborhood of (Gj) = (tEj) in {C{dft))'^. Since the measurements D = D, we 
have that L and L solve the same transport equation ([55]) while L = L = D/G on 
As /3 satisfies ((69)) . we deduce that L = L since the integral curves of /3 map 
any a; G fi to the boundary dil. More precisely, consider the flow 9x{t) associated 
to /?, i.e., the solution to 

e^{t)^P{e,{t)), e{Q) = xe^. (74) 

By the Picard-Lindelof theorem, (|74p admits a unique solution since /? is of class 
C^(f2). For a; e let x±{x) e 9fi and t±{x) > such that 

e^{t±{x)) ^ x±{x) £ on. (75) 

By the method of characteristics, the solution L to the transport equation ([55)1 is 
given by 

i(a;) =io(a;±(a;))e"^«*'^'^^^-^"^^''". (76) 

The solution L is given by the same formula since 9xit) — Ox{t). This implies 
Ej = Ej = Dj/L, j — 1,2. By the choice of illuminations, we have \Ej\ ^ 
due to (17^ and 7^ 0. Therefore, fc^n is uniquely determined by (15^ and thus 

tJ = CT. □ 



3.3 Vector fields and stability result 

Recall that B^if) is the flow associated with /3. From the equality 

- exit) = f mexis)) - p{0x{s))]ds, (77) 

Jo 

and using the Lipschitz continuity of /3 and Gronwall's lemma, we deduce the exis- 
tence of a constant G such that 

\Ox{t)-ex{t)\<Gt\\p-p\\co^^^, (78) 

when 9x{t) and 9x{t) are in Cl. The inequality (j78|) is uniform in t as all character- 
istics exit Cl in flnite time. 

To see higher order estimates, we deflne W := DxOx{t), which solves the equa- 
tion, W = Dx(3{6x)W, with W{0) = I. Define W similarly. By using Gronwall's 
lemma again, we deduce that 

\W-W\<Gt\\DxP-Dx^\\co^n) (79) 

when 9x{t) and dx{t) are in f2. Since /3 and {3 are of class G'^{Ct), then we obtain 
iteratively that 

\Di-'9x{t) - Di-'9x{t)\ < Gt\\Dxp - Dx^Wc^-i^n), (80) 

when 9x{t) and 9x{t) are in f2. 

Recall that fii is defined to be the subset of il by removing a neighborhood of 
each tangent point of dfl with respect to Co- 
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Lemma 3.7. Let fl be an open bounded and convex subset in with C"^ boundary. 
Let d > 2 and assume jS and P are C^ipl) vector fields which satisfy Ij69p . Restricting 
to rii, we have that 

- i+llc-Hs^i) + ^ hWc^-^ni) < - /3+||c<i-i(j2i), (81) 

where C is a constant depending on ^l. 

This lemma is similar to the lemma 3.8 in p^ and the lemma 4.1 in |4l, but uses 
a different proof. 

Proof. For a: G f^i, let 9x{t) and 9x{t) be two flows associated to vector fields (3 and 
respectively. Denote A := 9x{t+{x)) € and B := 9xii+{x)) G dfl. Without 
loss of generality, we assume i+(a;) < t+{x). We also denote C := 9a;{t-i-{x)) E H.. 
To simplify the writing, let S :— \\/3 — /3||c''(r2i)- 




Figure 1: Vector fields /? and (3 

We first want to show that the angle ZAxB is controlled by 

ZAxB < Ci6 + C2h, (82) 

for some Ci, C2. Indeed, by applying ([78| and sine theorem, we can see that ZAxC 
is bounded by CiS. Also notice that /3 satisfies ([M]) . Therefore, similar argument 
shows that, for any ^1,^2, the angle between the vector from x to 9x:{ti) and the 
vector from x to 9x{t2) is bounded by C2/1. Thus ZCxB < 6*2/1. This proves 

By the definition of fii, a neighborhood of the boundary point at which the 
tangent plane of dQ is parallel to <^o is removed, therefore, there exists a constant 
value (f)o > depending only on fli such that, for any 2^ G 0i > 0Oj where 
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(f)i is the angle between the vector xA and the tangle plane of dfl at A, as in Fig. 
[TJ Then by (|82p . when 6 and h are so small that 0q := 0o — CiS — C2/1 > 0, the 
extension of xB will intersect the tangent plane of dH. at A, with intersection point 
D. Then it is easy to check that 

ZABC > ZADx ^4)1- ZAxB >(t)'o>0. (83) 

The sine theorem gives that 

([781) directly implies \ AB\ = \x+ - S:+\ < C'S. 

Since /Sje C"'(17) and 917 is of class C"*, it is clear that ZABC and ZACS are 
functions with respect to x G il. By differentiating (|84l) and applying (f80| . we 

get higher order estimates 

||x+-£+||c<^-i(o,)<C"^. (85) 
To see the second part in (|8T|). we have that 

|Ci3| = m{d))ds = \p{e,{T))\{i+{x) - t+ix)), (86) 

for t+{x) < T < t^{x). Similar argument shows the estimate of i+ — 1+ in (fSTj) . □ 

Proposition 3.8. Lei > 1. Lei L and L be solutions to ()55|1 corresponding to 
coefficients (/3,7) and (/3,7), respectively, where (|69p /loZds /or &ot/i /3 and /3. 

Lei Lp = L|go and Lq = Lan, thus Lo,Lq g C"'(9r2). M^e also assume h is 
sufficiently small and is convex. Then there is a constant C such that restricting 

to rii 

\\L - L\\c'^-i(i-i^) <C'l|io||c''(90i)[||^ - ^11^-1(01) (87) 
+ II7 - iWc-^iUi)] + C\\Lo - Lollcd(aoi) 

The proposition and the proof follows proposition 4.2 in [4|. 

Proof. By the method of characteristics, L(x) is determined explicitly in (j76p . while 
L has a similar expression. We thus have 

\L{x)-L{x)\ <|(Lo(a;+(:r))-Lo(i+(x)))e-^o + *^'7(«.(^))'i^| (88) 
+ \U{x+{x)){e- ^o^'"'' _ g- 

Applying Lemma [3.71 '^e deduce that 

|L»f"^[Lo(a;+(x))-Lo(i+(a;))]| < ||Lo-Lo|| 0^-1(90) +C||io||c<'-i(an) 11/3-/311^-1(01) ■ 

(89) 

This proves Lo(x+(a;)) is stable. To consider the second term, by the Leibniz rule, 
it is sufficient to prove the stability result for /g*^"^^ ^{9x{s))ds. 
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Assume without loss of generality that < Then by applying (15(71) . 

we have 

Jo 

< C||7llc"(Oi)ll^ - i^Wcifii) + C||7 - 7llc"(f2i)- 

Derivatives of order d — 1 of the above expression are uniformly bounded since 
t+{x) £ C'^^^{ili), 7 has derivatives bounded on il and 9x{t) is stable as in 

It remains to handle the term 

v{x) := / j{0x{s))ds. 
Jt+ix) 

/3 and 7 are of class C"'(i7i), so is the function x — >■ ^{9x{s)). Derivatives of order 
d — 1 of v{x) involve terms of size t^{x) — t^{x) and terms of form 

D^(^i+Di-'-"^j{Ui+)) - t+Dt'-"^^{Ut+))y < m < d - 1. 

Since the function has d—1 derivatives that are Lipschitz continuous, we thus have 

\Di-'v{x)\<C\\t+-i+\\c--Hn,)- 

The rest of the proof follows Lemma [3771 □ 

Now we can prove the main stability theorem. 

Proof of Theorem \2.S\ The first part follows directly from ((55)) and Proposition 
13.81 This also provides a stability result for Ej = Dj/L as L is no n- vanishing. 
By choosing the boundary illuminations close to the boundary conditions of CGO 
solutions, (|7ip and ([7^ imply that Ej is non-vanishing since the CGO solutions are 
non- vanishing. Thus ([T5|) gives the stability control of fc^rt and thus a. □ 

3.4 Stability with 6 complex internal data 

Rather than applying the characteristics method to ([55]) . we can rewrite (f55|) into 
matrix form by introducing more internal measurements. We first construct proper 
CGO solutions. Let j = 1,2,3 in this section. We can choose unit vectors Co and 
tJq, such that Co ' Co = Co ' Vo — ^ ^^'^ {Col are linearly independent. Also, choose 
(Cf,C^') and irji,ri) such that |C| := |C^I = IC^I, 

C9 ■ ■ ■ 

" Vi^'^o^ ^^'^ ,i™ v^^vi- (90) 

We construct CGO solutions £'i,£^2 corresponding to {Ci,vi) and (C2:'72)- Let 
the boundary illuminations , G2 be chosen according to ([7T|) for e small enough. 
The measured internal data are then given by , Z?2 • Proposition 13.51 shows that 
the vector field defined by ([55]) satisfies that 
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\W'-L'Q\c^in)<j^^- (91) 

While Id is sufficiently large and L 7^ on f2, we obtain that the vector {f3^{x)} are 
linear independent at every x € f2. Thus matrix (/3^(x)) is invertible with inverse of 
class C"*(i7). By constructing vector- valued function r(x) G (C"'(il))^, the transport 
equation (1551) now becomes the matrix equation 



Vi + r(x)L^O. (92) 

Notice that r(a;) is stable under small perturbations in the data D := (_Dj,Z?2) ^ 
iC^m'', i.e., 

I|r - r||(c^-i(o))3 <C\\D- D\\^c^^n)r. (93) 



Assume is connected and Lq — L\qq is known. Choose a smooth curve from 
X G n to a point on the boundary. Restricting to the curve, (|92l) is a stable ordinary 
differential equation. Keep the curve fixed. Let L and L be solutions to ((92|) with 
respect to F and F, respectively. By solving the equation explicitly and (|93|). we 
find that 

\\L - L\\c^-iin) <C\\D- D\\^cHn)r- (94) 

Proof of Theorem HTM The first result in ([14]) is directly from The proof of 

the stability of a is exactly the same as in the proof of theorem 12.21 □ 

4 Acknowledgment 

The authors thank Professor Gunther Uhlmann for suggesting this problem and for 
helpful discussions. The authors also thank Professor Maarten de Hoop for helpful 
discussions on the inverse problem models. The work of both authors was partly 
supported by NSF. 

References 

[1] G. Bal and G. Uhlmann, Inverse diffusion theory of photo-acoustics, Inverse 
Problems, 26 (2010), 085010. 

[2] M. A. Biot, Theory of propagation of elastic waves in a fluid-saturated porous 
solid. I-Low-frequency range, Journal of the Acoustical Society of America, vol. 
28, pp. 168-178, 1956. 

[3] M. A. Biot, Theory of propagation of elastic waves in a fluid-saturated porous 
solid. II-High-frequency range. Journal of the Acoustical Society of America, 
voL 28, pp. 179-191, 1956. 

[4] J. Chen and Y. Yang, Quantitative photo- acoustic tomography with partial data, 
Inverse Problem, 28 (2012), 115014. 

[5] D. Colton and L. Paivarinta, The uniqueness of a solution to an inverse scatter- 
ing problem for electromagnetic waves, Arch. Rational Mech. Anal. 119 (1992), 
59-70. 



16 



[6] S. R. Pride, Governing equations for the coupled electro-magnetics and acous- 
tics of porous media, Phys. Rev. B, 50 (1994), 15678-15696. 

[7] S. R. Pride and M. W. Haartsen, Electroseismic wave properties, Journal of 
the Acoustical Society of America, vol. 100, no. 3 (1996), pp. 13011315. 

[8] A. Tliompson and G. Gist, Geophysical applications of electro-kinetic conver- 
sion, The Leading Edge, vol. 12 (1993), pp. 11691173 

[9] J. Sylvester and G. Uhlmann, A global uniqueness theorem for an inverse 
boundary value problem, Ann. of Math., 125(1) (1987), pp. 153-169.1 

[10] K. L. WilHams, An effective density fluid, m,odel for acoustic propagation in 
sediments derived from Biot theory, J. Acoust. Soc. Am. Volume 110, Issue 5, 
pp. 2276-2281 (2001) 

[11] Z. Zhu, M. W. Haartsen, and M. N. Toksoz, Experimental studies of electro- 
kinetic conversions in fluid-saturated bore-hole models. Geophysics, vol. 64, no. 
5, pp. 1349-1356, 1999. 

[12] Z. Zhu and M. N. Toksoz, Cross hole scismo-clcctric measurements in bore-hole 
models with fractures, Geophysics, vol. 68, no. 5, pp. 1519-1524, 2003. 

[13] Z. Zhu and M. N. Toksoz, Seismo-electric and seismo-magnetic measurements 
in fractured bore-hole models, Geophysics, vol. 70, no. 4, pp. F45-F51, 2005. 



17 



